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Abstract 



The vacuum two-point correlation function is calculated beyond the Gaussian 
approximation during the second order phase transition. It is found that the 
correlation function is dominated by the Gaussian term immediately after 
the phase transition but later is taken over by non-Gaussian terms as the 
spinodal instability continues. The non-Gaussian effects lead to larger size 
domains and may imply a smaller density of topological defects than that 
predicted by the Hartree-Fock approximation. 
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The phenomenon of continuous phase transitions is common to many diverse areas of 
physics. In particular phase transitions have been studied extensively in condensed matter. 
However the dynamical process of the onset of phase transition and temperature dependence 
or more generally its time-dependence is less well-understood. Topological defects may be 
formed during spontaneous symmetry breaking phase transitions. A theoretical understand- 
ing for defect formation in such phase transitions was provided by Kibble and Zurek 
on the basis of causality principle. 

Since then, it is argued that the formation of vortices may be observed experimentally 
in liquid 4 He || and 3 He ||. It is well known that liquid 4 He has a A-point phase transition 
to a condensed phase, superfluid 4 He, at T = 2.17 K. For liquid 3 He the phase transition to 
a superfluid state occurs at 2 mK. The transition in 3 He is more like the case of supercon- 
ductivity in metals where pairs of electrons are coupled to form loosely bound bosons in 1 S 
state, except that now 3 He is a fermion and the pairs are in 3 P state. However it appears 
that experiments do not give a definite result on the vortex formation. 

The early Universe might have undergone a sequence of spontaneous symmetry breaking 
phase transitions Heavy-ion collisions at relativistic energies lead to a quark-gluon 

plasma, which is believed to mimick the early stage of the Universe. In the process of 
cooling this plasma, the system will go through a chiral phase transition. It is likely that 
the phase transition may lead to the formation of disoriented chiral condensate (DCC) ||. 
The evidence of DCC would be emission of large number of coherent pions from the decay 
of such condensates. 

Current methods to study the dynamical formation of defects from the onset of phase 
transitions are classified largely into the time-dependent Landau-Ginzburg (TDLG) equation 
and quantum field theory (QFT) [j3-12]|. In the former the dynamics of phase transitions 
is described by the classical TDLG equation with a stochastic noise, whereas in the latter 
the dynamics is described by the quantum field equation with or without the noise. QFT 
is more fundamental than TDLG equation in the sense that QFT is the best tested, micro- 
scopic theory of nature. But QFT has been limited to, more or less, the Hartree-Fock or 
Gaussian approximation due to technical difficulty However, during the second order phase 
transition via a quench, the soft (long wavelength) modes grow exponentially. So it is likely 
that nonlinear terms would become comparable to the Gaussian contribution during the 
phase transition and spinodal instability quickly turns the dynamics of the problem to be 
highly nonlinear and nonperturbative. It may then be anticipated that large non-Gaussian 
fluctuations will drive the domain size. 

In this Letter we find the nonequilibrium quantum evolution beyond the Gaussian ap- 
proximation during the second order phase transition and calculate the correlation length 
from the two-point function. As a QFT model we shall focus on a real scalar field undergo- 
ing an instantaneous quench that describes the dynamical formation of domains. For that 
purpose we use the Fock basis obtained in Refs. |T2|-p^| to find perturbatively and system- 
atically the improved vacuum state during the quench of the second order phase transition. 
Our scheme is more systematic than an improved Gaussian state corrected by the second 
excited state |15] |, since we can find all the corrections to the Gaussian state to any desired 
order of the coupling constant using the standard perturbation method. 

The QFT model for the second order phase transition is described by the Hamiltonian 
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The quantum dynamics of nonequilibrium process is prescribed by the time-dependent func- 
tional Schrodinger equation 
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we rewrite the Hamiltonian ([[]) as 
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with a and /3 denoting {k, (±)}. In deriving Eq. (||) we neglect the odd power terms from the 
self-interaction whose expectation values vanish with respect to symmetric states because 
we are interested in improving the Gaussian state that preserves the parity. The mode- 
decomposed Hamiltonian ([]) consists of coupled quartic oscillators, the coupling among 
different modes arising from the last term. Because of the time-dependent mass squared 
the Hamiltonian (^) describes a system in nonequilibrium, the degree of nonequilibrium 
being determined by the rate of change of the mass. In order to handle this nonequilibrium 
system we use a recently introduced canonical method, the so-called Liouville-von Neumann 
approach dT^HH- 

The essential idea of Refs. []12-14}1 is to introduce the following time-dependent annihila- 
tion and creation operators for each mode that satisfy the quantum Liouville-von Neumann 
equation: 
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Then the Hamiltonian (^j) can be represented in terms of A a , A^ a and divided into 
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H(t) = H G (t) + (Xh 2 )H P (t), 
where H G is the quadratic (Gaussian) part and Hp is the quartic part: 
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Here we neglect purely c-number terms. The requirement that each A a and A\ satisfy the 
quantum Liouville-von Neumann equation 



ih^A a (t) + [A a (t),H G (t)]=0, 

ih^Al(t) + [Al(t),H G (t)]=Q, 



leads to the mean field equation for each mode 
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Further the standard commutation relations L4 Q (t), = 5 a ^ are guaranteed for all 

times by the Wronskian conditions 
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Then the Fock space for each mode consists of the Gaussian vacuum and excited states with 
quantum number n a 
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The Fock space basis satisfies the orthonormality condition G (m a , t\np, t) G = 5 a p5 mn . Up to 
a time- dependent phase factor each state \n a ,t) G satisfies the Schrodinger equation for the 
corresponding mode of H G . The wave functional for H G is a product of such a state for each 
mode. For instance, the exact vacuum state |0,£)g = I\ a |0 a , £)g f° r H G (t) is the Gaussian 
vacuum for the total Hamiltonian H(t) in the Hartree-Fock approximation. Other excited 
states of the Gaussian vacuum are the multiparticle states of modes, which are compactly 
denoted as 
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where {A/"} = (n\, • ■ ■ , n a , • ■ •) with all n a being nonnegative integers, A^^ 
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Quantum states beyond the Gaussian approximation can be obtained by treating Hp as 
a perturbation to H G . The perturbation .Hp excites and de-excites in pairs the multiparticle 
states, so we look for an exact quantum state of the form 



\{Af},t) = U[Al(t),A a (ty,t,\]\{Af},t) G . 
Then the Schrodinger equation (|3]) takes the form 



i^U(t, A) + [U(t, X),H G (t)] - XH P (t) 
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Using Eq. (|12|) and technically assuming that the time derivative does not act on A^ a (t) and 
A a (t), the operator U satisfies an interaction picture-like equation 



ih^U(t,X) = (Xh 2 )H P (t)U(t,X). 
The formal solution to Eq. (pi]) can be written as 



U (t, A) = Texp 



-Xh / H P (t)dt 
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where T denotes a time-ordered integral and A\(t) and A a (t) are treated as if they are 
constant operators. The operator U is a unitary operator, which is a consequence of the 
Hermitian H p and is easily shown from the formal solution (|20|). The unitarity of U preserves 
the orthonormality of ({Ai}, t\{Af}, t) = 5{m},{N}- The vacuum state is then given by 



\0,t) = U(t,X)\0,t) G = X: U m (t)\{Af},t) G . 
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Our vacuum state (21) is non-Gaussian in that it is a superposition of the Gaussian vacuum 



and its excited states flTJ], and further has kurtosis (higher moments) different from the 
Gaussian vacuum |U| . 

Before the quench (t < 0), the Gaussian vacuum for each mode is given by the solution 
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where Eq. (|T3| ) leads to the gap equation and yields the correct Gaussian vacuum energy 
|17|| . The contribution of non-Gaussian terms to the correlation function is of the order of 
(Xh) 2 /(2 ll Qf k ). Thus not only for the weak coupling limit (A mf) but also for the strong 
coupling limit A ~ m 2 , the non-Gaussian contribution to the two-point function is quite 
negligible. This justifies the validity of the Gaussian vacuum before the phase transition via 



5 



the quench. However, after the quench (t > 0), the soft modes (k 2 < m 2 ) begin to grow- 



exponentially due to m 2 (t) 
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The exponential growth continues until it reaches the 



spinodal line where (A^.)(^ k </5/ i k) from the self-interaction takes over the unstable 
Immediately after the quench we find approximately the soft mode solutions to Eq. (|H 
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where flf^ 



m. 



k 2 . These soft mode solutions after the quench continuously match 



those before the quench and become exact for the non- interacting case (A = 0). 

However, during spinodal instability the perturbation (Xh)Hp grows exponentially, be- 
comes comparable to the self-interaction term of Hq and will make a contribution to the 
correlation function as much as the Gaussian contribution. Since the soft mode solution 
(f23|), which dominates over hard (short wavelength) mode solutions, has an exponentially 
growing factor overall time-dependent factor, Hp(t) approximately commutes each other for 
two different times t' ^ t. Hence the formal solution (^) has a good approximation 



U[i]{t, A) = exp 
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Moreover, the higher order terms in an exponential form of the formal solution (see Appen- 
dices of Ref. [|TI|]) are suppressed by the expansion parameter Xh, which consolidates the 
validity of the operator (|24|) as far as soft modes are concerned. Using Eq. ([24]) and finding 
the most dominant contribution to the two-point correlation function up to order (Xh) 2 , we 
obtain 



(o,t|0(x,t)0(o,t)|o,t) = x;E- 

n=0 a 
d 3 k 



-ixny 



G (0,t\[j ^ P (t)rft,0 2 (t)] (n) |O,t) G e lk - x 

(Xh) 2 r d 3 k h ' 



ik-x 



+ 



(2tt) 3 2 J (2tt) 3 12 



ikoc , (25) 



where [0,B]r n -\ = [O, [O, J3]( n _n] with [O, B] = B. Note that all the terms linear in (Xh) 
vanish and at order (Xh) 2 the mode-mixing terms contribute quadratic terms of ip-k or to 
the scale dependent part, which are negligible compared to the quartic terms in Eq. (f25p. 

Finally, after inserting the solution (^) into Eq. and explicitly doing the integral, 
we obtain the two-point correlation function of the form 
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where G(0, t)'s are slowly varying functions of t. The leading and sub-leading terms of the 
second integral of Eq. (p5|) cancel each other and the third-leading term is the second term 
of Eq. (0). The first and the second terms of Eq. (|26| ) are the Gaussian and the non- 
Gaussian contribution, respectively, to the correlation function. Each term oscillates and is 
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bounded by an envelope determined by exponentially decreasing functions. The first term, 
the Gaussian term, has a larger initial amplitude but decays more rapidly than the second 
term, the non-Gaussian term. Therefore, immediately after the phase transition, the first 
term determines the domain size, which obeys the well-known Calm-Allen scaling relation 

II 




= \hr- ( 27 ) 



However, as the phase transition continues, the second term begins to dominate over the 
first term and leads to the non-Gaussian scaling relation 



£ng(*) = J— = v^g. (28) 

V m f 

A few comments are in order. First, as the time for spinodal instability goes on, higher 
order non-Gaussian terms begin to grow comparable to and eventually become larger than 
the Gaussian and lower order terms. A series of transitions of scale relation may be expected 
from the Gaussian one to the non-Gaussian ones of higher orders. In fact it cannot happen 
in a realistic scenario because of the sampling of the mean field around the true vacuum 
and the de-excitation of excited states at order (Xh) 2 or higher. The sampling of the mean 
field <fk(t) around the true vacuum after crossing the spinodal line stops exponential growth 
and limits the duration of spinodal instability. Also there will be equally excitation and 
de-exitation, thus making virtually impossible for a system to achieve a state with inverted 
population, the very high excitation having the largest population. In reality there will be 
an increased occupation of excited states near the Gaussian state and the very highly excited 
states will have a small occupation probability. Second, though we focus on the real scalar 
field, the result above can be readily extended to a scalar field with O(N) symmetry before 
the phase transition. The two-point function of a complex scalar field is essentially the same 
as that of the real scalar field |12|]. We thus anticipate that other topological defects, such as 
monopoles, strings, and vortices, formed from the scalar field with the appropriate internal 
O(N) symmetry may have the same correlation length as domains. The non-Gaussian effects 
thus may decrease the defect density. 

Finally we conclude with the physical implication of non-Gaussian effects on domain 
growth, DCC and the density of topological defects. The RHIC at Brookhaven and the 
LHC at CERN in the near future will produce heavy ion beams at very high energies that 
will lead to deconfinement of quarks present inside individual nucleons. It is anticipated that 
a quark-gluon plasma will be formed and approach equilibrium over a time ~ 10~ 23 s and 
then will evolve rapidly far from equilibrium through a chiral phase transition. It has been 
suggested that the chiral phase transition will be accompanied by the formation of DCC ||. 
From a theoretical point of view it was argued that the linear a model may be a reasonable 
model, that includes some of the symmetry properties such as chiral symmetry of QCD, to 
establish the extent of DCC. Rajagopal and Wilczek have advocated that the exponential 
growth of soft modes of the linear a model may lead to a large correlation and hence a large 
domain of condensate [0 . This mechanism has similarities to the spinodal decomposition in 
QFT. The numerical simulation based on the Hartree-Fock approximation, however, results 
in smaller (1 ~ 2 fm) DCC regions than the anticipated (> 10 fm) DCC regions [19] . 
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On the other hand, we find that the non-Gaussian effects increase the domain volume by 
an order of magnitude. This implies that a chiral phase transition may result in the formation 
of a large DCC domain, the decay of which eventually leads to a profuse production of 
coherent pions. However any extraordinary observation should be investigated carefully to 
find if there is a linkage to the DCC. The present study of a real scalar field serves only 
as a simple model that suggests the role of non-Gaussian effects on domains, and implies 
that the density of topological defects may be small. A realistic study of DCC formation 
during the chiral phase transition would require working with the linear a model. It may be 
anticipated that long wavelength modes will play a similar role to enhance not only the size 
of domains in quark-gluon plasma, condensed matter and the early Universe but also the 
density perturbation from the second order phase transitions in inflation scenarios. These 
studies are in progress. 
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